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Abstract 

A new method of constructing an infinite set of exact eigenfunctions 
of 1-d interacting spinless Fermionic systems is derived. Creation and 
annihilation operators for the interacting system are found and thereby 
the many-body Hamiltonian is diagonalised. The formalism is applied to 
several examples. One instance is the theory of Jack polynomials. For the 
Calogero-Moser-Sutherland Hamiltonian a direct proof of the correctness 
of the asymptotic Bethe Ansatz is given. 



1 Introduction 

The study of one-dimensional integrable models of interacting particles has a 
long history going back to Bethe |Bet31| . In physics there has been a renewed 
interest in one-dimensional integrable systems in the last years in the study of 
cold atom gases and Bose-Einstein condensates |FRZ04[ ICZ04a[ ICZ04b| . 

The standard way of constructing eigenfunctions is Bethe's Ansatz. The 
crucial condition for Bethe's Ansatz to be successful is that the two-body scat- 
tering matrix Sabikj — kk) of a particle of type a with momentum kj with a 
particle of type b with momentum kk fulfills the Yang-Baxter (star triangle) 
equation, which serves as the starting point for the algebraic Bethe Ansatz. 

Imposing periodic boundary conditions for the N particle wave function, 
from Bethe's Ansatz one (or a set) of Fredholm integral equations for the density 
of states is obtained. They are referred to as Bethe Ansatz equations. 

Applied either to 1-d quantum mechanical, or to 2-d classical lattice the- 
ories, Bethe's Ansatz has been extraordinarily successful. For one reason, be- 
cause many of the most important lattice models, as for instance the Hubbard 
model, Heisenberg model, Ising model etc. only nearest neighbor interaction is 
assumed. Applied to continuous models Bethe's Ansatz is in its simplest form 



*heinerich.kohler@uni-due.de 



1 



constrained to particles with (5-interaction |LL63[ lYanSTi IGau66| . being the only 
strictly local interaction. 

The problems, connected with non-local interactions in continuum models, 
were partly overcome by the asymptotic Bethe Ansatz (ABA) . It was introduced 
by Sutherland [SutTlai ISut71b| in order to obtain thermodynamical quantities 
for Calogero-Moser-Sutherland (CMS) models. 

The basis assumption of ABA is that the Bethe Ansatz equations still hold 
for non-local interactions as long as the iV-body ^-matrix factorizes into a 
product of 2-body S'-matrices. Using the ABA hypothesis it is therefore possible 
to obtain thermodynamical quantities without detailed knowledge of the wave 
function in the interacting region. The correctness of this assumption has been 
proven for the trigonometric CMS model, where a complete set of eigenfunctions 
can be constructed with Jack polynomials [Mac95[ IFor92[ IFor931 IFor95| . In 
that specific model even some thermodynamical correlation functions could be 
calculated |Ha95| . 

Although the correctness of ABA has been proven also in some other cases 
|Kaw92] , the different treatment of local and non-local interactions within the 
Bethe Ansatz is unsatisfactory. Moreover ABA yields no clue of how to construct 
eigenfunctions of systems with non-local interaction. Therefore in this work I 
wish to put forward an approach, alternative to Bethe's Ansatz which treats 
local and non-local interactions on the same basis. 

The exact A'^-body wave function is not constructed via an Ansatz and by 
adjusting parameters but by the successive application of a creation operator 
onto the vacuum ground state. Thereby an integral representation for an ar- 
bitrary TV particle eigenfunction is obtained. The explicit construction of this 
creation operator and its corresponding annihilation operator for certain inter- 
action potentials is the main result of this work. 

The basic idea roots on the following observation: a class of multidimensional 
integration formulae - some of them are well known for a long time - allow for 
a natural interpretation as Fermionic creation operators of a one-dimensional 
many-body Hamiltonian. These integrals share the common property that the 
integration domain /j^j^-' of a set x' of N integration variables , 1 < n < A'^ is 
defined by the interlacing condition 

Xi > x[ > X2 > ■ ■ ■ > XN > x'j^f > Xj^i^i . (1) 

The integral itself is therefore in general a function of a set x of iV+ 1 arguments 
Xi. Due to condition ([1]) this function vanishes, whenever two arguments are 
equal. Thus by construction it obtains the nature of a Fermionic wave function, 
when the A -f 1 arguments are interpreted as particle positions. One instance 
is the following version of the Dixon-Anderson integral |Dix05[ IAnd91| 

A , M Af W N N+1 

^a(x,x') = ;^5^nn(-^--;-)^n n (-;--ov(3) 

^N+l l^J i=l 1 = 1 1 = 1 i=j+l 
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It might be considered a representation of a power of Vandcrmonde's determi- 
nant 

N+l 

AAr+i(x) = (x„ - Xm) (4) 

n<m 

as an integral over Ajv(x'), if both sides are multipHed with A^^j'^(x). In 
Ref. |GK02] group integrals of the form 

</.^^(k,x) = / dfi{U) exp{-itYU-^xUk) , (5) 

Jueu{N) 

were considered. Here U{N) is a classical Lie group manifold and x, k are 
diagonal N x N matrices, and dii{U) is the Haar measure of the group. The 
parameter A depends on the group under consideratiorQ. It was found that 
0^'' (k, x) can be constructed recursively in N as 

•^An i=l 

where /eat+i is the additional element of the new {N + 1) x {N + 1) matrix k+ 
on the left hand side. Okounkov and Olshanski |0097) found an integration 
formula for symmetric Jack polynomials J^"^ (n, x) in N variable^ 

j(,'|i(n,x) cx l^^l[dx',f,xMJ^^^\ny), (7) 

which relates Jack polynomials with + 1 arguments to Jack polynomials with 
N arguments to the same partition n = {ni, . . . , njv}, rii e N. 

Eq. ([2]), Eq. ([6]) and Eq. ([7]) share the same structure: An integral of a func- 
tion f oi N variables multiplied with an integration kernel /iA(x, x') reproduces 
the same function / with the number of variables increased by one. I will embed 
these instances in a general framework and show that the proper generalization 
of /za(x, x') has a most natural interpretation as the coordinate representation 
of a particle creation operator . 

Whereas some of the integral formulae arising from the creation operator 
are already known, the construction of the annihilation operator is completely 
new. The resulting integral representation of the annihilation operator yields 
interesting new integral identities, which might be useful in applications. 

This paper is the first in a series, dealing with the construction of creation 
and annihilation operators in 1-d integrable systems. It focuses on spinless 
Fermions. Bosons and spin 1/2 Fermions will be addressed in separate publica- 
tions. 

The paper is organized as follows. The first section supplies a precise de- 
scription of the problem, whose solution will be presented in Theorem [TJ A 

is related to the parameter (3 of IGK02I by A = /3/2 - 1. 
^ A is related to the parameter a of MacDonalds book byA = l/o — 1. 
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discussion of Theorem [T] follows. In Sec. [3] the results of Theorem [T] are illus- 
trated in several applications. Proofs will be given in Sec. H] 



2 Statement of the result 

We consider a Hamiltonian, describing one dimensional non-relativistic spinless 
interacting particles with mass 1/2 

dxdx'iP^ (a;)V'^ {x')V{x - x')i}{x')il}{x) . (8) 



The integration domain U, is the real axis or a compact interval. For the first 
choice the spectrum of H will be continuous. This has the implication that the 
eigenvalue problem for H might have no solutions in the Hilbert space of 
square integrable -functions in M. In this case we look for solutions in an 
enhanced (rigged) Hilbert space, which allows for eigenstates of H which are 
not normalisable. In other words, we consider the Gelfand triplet V C C V , 
where V is the space of test functions and V is the space of distributions, dual 
to v. In particular V includes eigenfunctions, which behave in the asymptotic 
limit as plane waves (scattering solutions of the Schrodinger equation) . For this 
behavior in the asymptotic limit of the wave function we introduce the term 
scattering boundary condition (SBC). We consider 

^ _ f M , for scattering boundary conditions (SBC) , 

[[0, L] , for periodic boundary conditions (PBC). 

If we denote the one-particle Hilbert space by C (il), we have 

C (il) / ' ^'^'^ scattering boundary conditions (SBC) , ,^qs 

, for periodic boundary conditions (PBC). 

In Eq. ([8]) i/^ix) and ip'' (x) are Fermionic creation (annihilation) operators obey- 
ing the anticommutation relation 

{iPix),'4>'<ix')} = V(a;)V'^(a;') = Six - x') . (11) 

For a fixed number of particles N, the Hamiltonian can be written in first 
quantization in coordinate representation as 

ij^(x) = - E + E ^(^^ - ' (12) 

where x denotes the set of particle positions x = {xi, . . . , xn}- The Hamiltonian 
is the direct sum H = 0^^o acting on the Fock space H — U'j^^q'Hn, where 

'Hn — {V''(x) G ^ (^^) IV'(^) completely antisymmetric} , (13) 
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and £ {^^) = {^)- On Hn a scalar product is defined by the A^-fold integral 

(V-ivl^w) = / d^[x]V^(x)^jv(x) , (14) 

where the infinitesimal volume element [x] — Il^=i introduced. 

Our goal is to map the Hamiltonian as given in Eq. ([5]) onto the quadratic 
form 

H = J2e{k)alak , (15) 

k 

where Uk (aj^) annihilates (creates) a particle with quasimomentum k. The 
Hamiltonian acts on the complete Fock space and is block diagonal in the 
basis of eigenstates of the particle number operator 

In contrast, the creation and annihilation operators define a map which does 
not conserve particle number 

afc : TCn Ti-N-i , afclV'w) € TCn-i , 

al : Hn Hn+i , o^lV'Jv) e Hn+i ■ (17) 

This mapping reads in configuration space 

V'iv+i(x) = / d^[x']4(x,x')^jv(x') , 

^jv-i(x) = / d^[x>fc(x,x')^jv(x') . (18) 

This means, in configuration space the creation operator is an integral opera- 
tor whose kernel is a complex valued function of two sets of coordinates x = 
{xi, . . . , xn+i} and x' = {x'l, . . . , x'j^}. We call a|(x, x') creation function. The 
annihilation operator is an integral operator whose kernel afe(x, x') is a complex 
valued function of x = {xi, . . . ,a;jv-i} and x' = {x[, . . . ,x'j^}. We call afe(x, x') 
annihilation function. 

The creation and annihilation operators are defined by the basis independent 
commutator relation 

[H,ak] = -e(k)ak 

[iJ,4] = e{k)al. (19) 
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They read in coordinate representation 

= [i7w_i(x)+e(fc)] / d^[x>fe(x,x')Viv(x') 



= 



d^[x>fc(x,x')i/Ar(x')VAr(x') , 

[i/w+i(x)-6(fc)] / d^[x']4(x,x')^jv(x')- 

/ d^[x']4(x,x')ifAr(x')^Ar(x') . (20) 

This translates into a set of partial differential equations for the creation func- 
tions a|(x, x') and the annihilation functions aA;(x, x') 

[ifAr_i(x) - i7jv(x')] afe(x,x') = -e(fc)afc(x,x') , 
[i/w+i(x)-ffw(x')]4(x,x') - e(fc)4(x,x'). (21) 

With a set of operators fulfilling the commutator relations ([TOl) simultaneous 
eigenfunctions of the Hamiltonians Eq. p^)) and Eq. to the eigenvalue 

£' — ^(fcn) can be constructed by 



IViv(k)) = n4jO) . (22) 

n=l 

Since two operators which have the same eigenfunctions to the same eigenvalue 
are equal, the conditions specified in Eqs. (fTO|) to (|2ip are sufficient to prove 
that the Hamiltonians in the both forms as given in Eq. ([5]) and in Eq. (fT5|) are 
equal up to a basis rotation. 

Restricting ourselves to Fermions, we require the wave function (x|';/'7v(k)) € 
T^AT to be completely antisymmetric in two sets of arguments k = {fci, . . . , fcAr} 
and X = {xi, . . . , xtv}, thus the quasiparticle creation and annihilation operators 
have to obey the Fermionic anticommutation rules 

r \x f^C"- ' SBC , 

= I PBC. (2^) 

We introduce the following notation: ^„(k, x) denotes always a function 
with n € N arguments {fci, . . . , kn} and with n arguments {xi, . . . , x„}. 

We now state the main result of this section as a theorem. It states a suffi- 
cient condition on the interaction potential V{x) for the existence of the above 
defined creation and annihilation operators and gives an explicit construction 
of the eigenstates. 

Theorem 1 (Spinless Fermions) 1. Let f{x) be an antisymmetric func- 
tion satisfying the condition 

f{x)f{v) + f{x)f{z) + f{v)f{z) = const. , 

x + y + z = (24) 
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andF{x) = dx' f{x'), such that F{x) = F{—x), then annihilation func- 
tions afc(x, x') and creation functions a[.(x, x') satisfying the conditions 
specified in Eq. \21\) are given by 



ai;(x, x') = exp 



Af-l 



4(x,x') 



'i<m n,m 

/N-1 N \ 

\n—l m—1 / 

) + Y,F{x,,~x'^) 



exp 



N+l 



N /N+l N \ 

n<m \n— 1 m—l / 



(25) 



^. T/ie interaction potential in the Hamiltonian is related to the function 
f by 

V{x) = f^{x) - f{x) +const. (26) 



3. With the complex functions afc(x, x') and a|,(x, x') defined in Eq. 

eigenfunctions V'jVitiC*, x) to the N± 1 particle Hamiltonian Eq. can 
he constructed recursively as 

V'Ar+i(k,x) = (x|4^^JV'A'(k)) 



CAr(k) 



d^[x']4^^^(x,x')V^(k,x') , (27) 



V'Ar-l(k,x) (x|afc„|V'Ar(k)) 

= RVNCN-i{k) f d^[x>fe„(x,x')^Ar(k,x') ,(28) 

where R — (27r)~^ for SBC and R = L^^ for PBC. The integration do- 
mains //jf^'' and /q^^ are defined by 



/ d^[x'](...) = / dx', / dx'^... dx'^{...) 



(29) 



fXN-l 



p poo pXi pXN-1 

/ d^[x'](...) = / dx'J dx',... rfxW (...), (SBC) 



d-[x'](...) 



XI <J X2 



dx'^{...) , (PBC) 



The normalization constant CAr(k) is coordinate independent. For the 
most important potentials the explicit value will be given below (see Propo- 
sitionWil. 
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4- The dispersion relation of the quasimomenta is quadratic: 

e(fc) = fc2 . (30) 

The functions {21^ and \2l^) are eigenf unctions to the N± \ particle Hamil- 
tonians (I12p and to the center of mass momentum operator 

P^±i(x) ^ — , (31) 

n— 1 ^ 



such that 

i?jv±i(x)V'JV±i(k,x) = [J2kl]^N±iik,x), (32) 



/N±l \ 

PAr±i(x)V'jv±i(k,x) = fcri -0Af±l(k,x) 



(33) 



5. FFif/i respect to the scalar product Eq. i/ie orthogonality relation 

det[27r(5(fc^ - %)]i<jj<Ar , SBC 

l<ij<A' 



(^.(k')i^.(k)) cx { det[L4,..]::: ; ' pbc. (^4) 



The complete proof of Theorem [T] will be given in Sec.[31 Here we discuss some 
of its consequences. 

N-particle wave function As a corollary to Theorem [1] every iV-particle 
wave function can be written as a N{N — l)/2-fold integral as follows 



x("))) . (35) 



The coordinate sets x and x(^+^) are identified. Eq. ([35]) is obtained by iterating 
Eq. dlZl). 



Condition on the potential The functional equation is a special case 
of the functional equation 

/(a^)/(y)+/(a;)/W+/(y)/(^) = W(x) + W{y) + W{z) 

x + y + z ^ , (36) 

which was found by Sutherland [Sut71a[ ISut75j to be the necessary condition 
for a product wave function 

N 

^/'iv(0,x) = ]Jexp[P(x, -x,)] (37) 
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to be an eigenfunction to eigenvalue zero of an A'^-body Hamiltonian with two- 
body interaction only. The general solution of ((36)) was found by Calogero 
|Cal75bl[Cal75al . It is given by 

fix) = \ax\u,Lu') + l3x , (38) 

where ({x\uj,u!') is the Weierstrass ^-function |AS72| and uj,uj' are the two peri- 
ods of the corresponding Weierstrass elliptic function P. The right hand side of 
is determined by the Frobenius-Stickelberger equation for the Weierstrass 
C-function [FS80] 

iax)+ay)+az)f = -cix)-c{y)-c{z) 

x + y + z = , (39) 

to 

W{x) = -l{Xf'ix) + fix))-^. (40) 

The condition 3W{x) = const, therefore requires 

-^(A/'(x) + /2(x)-A/3) = const. (41) 

as well. Eq. (|^T|) can be integrated. The most general simultaneous solution to 
Eq. (gl]) and Eq. ^ is, with const. = -z^A^ 

f{x) = zAcoth(z(a; + a;o)) , (42) 

for arbitrary complex parameters xq, z and A. Requirement of a real potential 
restricts the values of z = a + ifc to be either real z — a oi purely imaginary 
z — ib and A to be real. These solutions are obtained from Eq. ([38]) by setting 
(3 — > a^A/3 and uj' — > i7r/2a, a; — > cxd or w — > 7r/26, uj' — > oo. A series of 
potentials can be derived by taking various limits. The most important cases 
are listed in Table [1] They are the interactions of the trigonometric (I), rational 
(II) and hyperbolic (III) Calogero-Sutherland-Moser (CMS) type. Remarkably 
also the sign-function respectively the (5-distribution are obtained from the 
shifted hyperbolic CMS Hamiltonian in the limit a — )■ oo, A — > with aA = c 
finite. 



Structure of creation and annihilation functions The creation (annihi- 
lation) functions as defined in Eq. ([25]) are strictly speaking also functions of 
the particle number N. We suppressed this obvious N dependence in order 
to unburden notation. However, to describe the relation between Uk and 

it is useful to indicate the N dependence by the symbols a^,^'' (x, x') and by 
al!'^\'x,x.'). Then we have 



a 



(x',x) = ar>M . (43) 
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/(^) 


V{x) 


F{x) 


const. 




(I) 


Xbcot{bx) 


6"A(A + 1) 
sin^(a;6) 


A In sin(6x) 


b^X^ 


limb_^o (I): 


(11) 


A 

X 


A(A + 1)^ 


Alnlxl 





b^ia (I): 


(III) 


aX coth(ax) 


a^A(A + 1) 
sinh^ (xa) 


A In 1 sinh(aa::) | 


-a2A2 


x^x + iTT/2 (III): 


(IV) 


aX tanh(aa;) 


a'XiX + l) 
cosh (a;a) 


A In 1 cosh(aa;)| 


-a^A^ 


lima^oo (IV): 


(V) 


c sgn(x) 


-2c5{x) 


c\x\ 





Table 1: List of interaction potentials for which creation and annihilation oper- 
ators can be constructed. They are based on the addition theorem for the cot 
function: cot(x)cot(y) = cot(x + y)(cot(x) + cot{y)) +1. 



The general structure of the creation (annihilation) function aj,(x, x') (afc(x, x')) 
factorizes into a k independent part aj(x, x') (ao(x, x')) and a k dependent part. 
The k dependent part is for all types of potentials a product of plane waves of 
the primed and the unprimed coordinates. They may be considered arbitrary 
eigenfunctions of the non-interaction V(x) = many-body Hamiltonian. On 
the other hand, as mentioned above, the basic ingredient of the fc~independent 
part is the product function tp]y{0,K) = JJ-^ - exp[F{xi — Xj)], which is an 
eigenfunction of the interacting system to eigenvalue zero. Using this general 
structure one might extend Theorem[l]to systems, where translation invariance 
is broken by an external potential, but whose ground state wave function can 
yet be written in the product form of Eq. (|37p . 



Boundaries of the integral operator Since we wish to create Fermionic 
eigenstates we cannot take ft^ as integration domain in Eqs. ([?7|) and ([^5]) . 
Because the integrand is by construction completely antisymmetric in the N 
integration variables x' this would yield the integral unavoidably zero. By the 
same argument the integral becomes zero, whenever two integration variables 
have the same integration domain. The integration domains //if^' and ij^^^ as 
chosen in Eq. (|29|) guarantee that the wave functions '0iv±i(k, x) are completely 
antisymmetric. By the definition of //jf and ij^^^ no ordering of particles is 
introduced. The newly created wave functions are defined in any of the (N ± 1)! 
disjunct regions of fl^^^. 

It is non-trivial to see that by the introduction of the domains I^^^ and ij^^^ 
the commutator relations in coordinate representation ([20]) are not affected. 
Since Hn± (x) acts also on the boundaries it has to be shown explicitly that no 
additional boundary terms occur. For periodic boundary conditions this yields 
a restrictive condition onto the form of the creation (annihilation) functions. 

Proposition 1 (Periodicity) For PBC the creation (annihilation) function 
a[.(x) ('afe(x) ) itself must be periodic with period L. 
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The proof of this statement is given in App.lXl Proposition[T]essentially fixes the 
boundary condition for aU potentials in Table [TJ Potential (I) is periodic with 
period -k /h. Therefore the asymptotic regime is never reached and only periodic 
boundary conditions L = ii/h are allowed for this potential. By Prop. [T] for 
potentials of type (II) to (IV) only scattering boundary conditions (SBC) are 
allowe41. Only for particles with ^-interaction both boundary conditions can 
be imposed. 



Analyticity at the Boundaries For the integrals ((27)) and (|28p to exist the 
creation and annihilation functions can have at most an integrable singularity 
in their integration domain. For potentials (I) to (III) this yields a restriction 
for the coupling parameter A to the range A G [— l,oo). This restriction is 
consistent with the well known fact that for CMS-Hamiltonians the minimal 
value of the coupling constant g — A(A + 1) is g = —1/4 |OP83| . For these 
potentials the behavior of -0^ (k, x) when two particles come close to each other 
can also be extracted directly from the creation and annihilation functions. The 
wave function vanishes with a typical power 

lim 0^(k, x) cx {x, ~ x,+i)^+i + 0{[x, - .T,+i)^+2) . (44) 
This shows again that A > — 1 must be imposed. 



Normalisation For potentials (II)-(V) with SBC the wave function is not 
normalisable. The constant CAr(k) can be evaluated by requiring that 'i/;7v(x, k) 
becomes the free particle solution in the asymptotic regime, i. e. in the regime 
where all distances Xi — xj, 1 < i,j < N are large. For potential (I) the wave- 
function is normalisable and the calculation is more involved. The result for all 
cases can be summarized as follows. 

Proposition 2 Depending on the potentials listed in Tahle\^the k-dependent 
normalisation constant CAr(k) is given by 



Cn(M) 



L 

,-Af(A+l 



-N N 



n 



{2i)^sgii{ki - k 



N+l) 



,A + 1 



r^(A + 1) 



i\ B {^\k, - kr,+i\ 

N 

kN+i\^ ih - kN+i) 



1=1 

N 



for (I) , 
for (II) , (45) 
for (III)-(V) , 



1=1 

where B{x, y) is Euler's beta-function. The proof of Prop.[2]is given in Sec 



Continuing potential (11) periodically leads to potential (I). 
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3 Applications 



The N{N — l)/2 integrals in the representation Eq. p5)) of the exact N par- 
ticle state can in some cases be evaluated exactly. However, many important 
properties of the wave function can actually be extracted from Eq. p5p without 
solving the integral. In the following we discuss free Fermions as well as the 
interaction potentials (I), (II) and (V) of Tab.[TJ 



3.1 Free particles 



Since it is instructive to see, how Theorem [T] works in the simplest case, we 
illustrate it first for free Fermions. A general iV-|- 1 particle state is given as the 
Slater determinant 



(x|^jv+i(k)) 



1 



= det[e^'="-]^<„„^<^_^^ ^ <^,(k,x).(46) 



This result can be derived from Theorem [T] by induction as follows: Assume the 
N particle states have been constructed, then we find the iV + 1 particle states 

by 



(x|V'^+i(k)) = (xi4„^jv^(k)) 



^''^^^ I d^[x']4 (x,x')det 



l<n,m<N 



C^iv(k) 



d^[x'l det 



l<n,m<N 



(N+1 N \ 

i/cjv+i ikN+1 ^ x[,r ■ (47) 

n— 1 ni — 1 / 



With the integration domain I^^'' given by Eq. (|29p the integral can be per- 
formed yielding again a determinant 



(x|^jv+i(k)) 



det [e''^""'"]^ 



<n,?Ti<JV+l ' 



(48) 

which is the desired result (compare with Prop. [5]). The action of on an N 
particle state is given by 

(xlVAr-i(k)) = (x|afc|V'w(k)) 



RC N- 


-i(k) 




-1)! 


RC 


-i(k) 




-1)! 



d^[x']afe(x,x')det 



l<n,m<N 



l<n,m<N 



N-1 



N 



exp j ik Xn — ik x'^, 



(49) 



m— 1 
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If we choose SBC we have to equip the quasimomentum k with a positive imag- 
inary increment ie for the x'l integration and with a negative increment for the 
x'pj integration. Then the integrals are convergent and yield again a determinant 



(x|V>Ar_i(k)) = 



1 Cjv-i(k) 



2n ^{N-l)\ 
2nS{k - fcl) 



2TrS{k-kN) 



iki — ik 



ikjsi — ik 



iki — ik 



ik 



N 



ik 



(50) 



From Eq. 



50)1 it follows 

afclV'Ar(k)) 
(xlafclV'Ar(k)) 



, if fc 7^ fci, i = 1 . . , 
Cjv-i(k) !_ 



N 



ik) V(^-l)! 



det 



l<n< W-1 



for A: — kn 



(51) 



which is the desired result (compare with Prop. [2] with j ~ N). 

For PBC we require vjv^C^, x) to be a periodic function with period L. 
Therefore k = x'^' '^here n is a set of N integers G Z. For the action of 
the creation operator nothing changes as compared to SBC. For the action of 
the annihilation operator we find that the (5-distribution is substituted by a 
Kronecker-delta 

(5(fc,-fc,) ^ ^K,k, ■ (52) 
This completes the construction for free spinless Fermions. 



3.2 Particles with 5— interaction 

For spinless Fermions the J-interaction is invisible and the wave function be- 
comes identical with the wave function of free Fermions as given in Eq. (|46|) . 
This result is quickly derived using Theorem[TJ Using the free solution ■0/^'' (x, k) 
for ^jv(x, k) in Eq. ([77|) . we obtain 

(x|V'^+i(k)) = 



(xi4ivA,(k)) 

CAr(k) 



d[x']4(x,x')det 



Kr 



.(53) 
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For the creation function a| (x, x') wc use Eq. (|25p and extract F{x) — 
from Tab. [1] 



(N+l N \ 

71—1 m— 1 / 



(54) 



exp 



n>m 



It is readily seen that the exponent in the second hne of Eq. ((54)) drops out 
completely for an arbitrary order of the particles. Eq. (|53p reduces to the free 
particle expression Eq. (|T7l) . The corresponding result is obtained for the action 
of the annihilation operators. 

3.3 Trigonometric Calogero—Moser— Sutherland system 

For the interaction potential (I) the model is called trigonometric CMS-model. 
As was pointed out in Sec. [2] the Hamiltonian has to be considered with PBC. 
First we recall some facts of the model. 

The eigenfunctions of the trigonometric CMS-models can be written as 

N 

V'Ar(n,z) = Vw(0,z)]^zf J^^'^(n,z) , 



1=1 

N 



■0^(0, z) 



aA+1 



-^(w-i) 



(55) 



i=l 



where the arguments Zi are related to the particle positions Xi by the relation 

Zi ^ exp (2TTixi/ L) . (56) 

Here -0Ar(O,z) is the ground state wave function of Hn and J^^^(n,z) is a 
symmetric polynomial in N variables Zi labeled by a partition n — {ni, . . . , un) 
of integers ni < n2 < . . .un- These polynomials are called Jack polynomials 
and were extensively studied |Mac95[ ISta89j . The parameter A is related to 
the parameter a of McDonalds book [Mac95| by A = 1/a — 1. The center 
of mass momentum i^T is a real parameter. Jack polynomials are defined as 
eigenfunctions of the operator H'j^ ^, which is obtained by adjunction of iJjv 

with the ground state wave function "0^(01 z) times the Galilean boost Hi^i 

-1 / 



N 



^jv(0,z)]^; 



K 



1=1 
N 



( 



N 



H 



N 



0jv(O,z) J|; 



E 

,2=1 

+2(A + 1)^ 



_d_ 

' dzi 

N 



+ [(\ + l)(N -l) + 2K]Y,. 



i<j 



d 



d 



Zi — Zj \ dzi dzj 



(57) 
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The ground state energy E^^'^^ is given by 

i;^''^^ = ^{N + l)N{N + + NK{K +{\ + l){N -\)) . (58) 

For K = and A = this is identical to the ground state energy of free spinless 
Fermions. The wave function ip^i^, z) in Eq. (|55p vanishes with the power A + 1, 
when two particles come close to each other, but it is not antisymmetric under 
interchange of two particles. Rather it obtains a phase 7r(A+ 1) under the action 
of the permutation operator P 

Primlpixi, . . . ,Xn, ■ ■ ■ ,Xm, ■ ■ ■ ,Xn) = 1p{xi, . . . , Xm, ■ ■ ■ , Xn, ■ ■ ■ , Xn) , Vn, m . 

For A even '0Ar(n, z) is Fermionic, for A odd it is Bosonic. For arbitrary 
real A it is a wave function with anyonic statistics |Ha95| . A Fermionic wave 
function can always be obtained from Eq. (j55p by the substitution 



N N 

K^'i^m^-"^'""' - |A^(z)|^A.(z)n^r'""^^/^ . (60) 

i=l i=l 

Jack polynomials depend strictly speaking also on the center of mass momentum 
K. In the following we fix the Galilean boost to K = —^^^{N— 1) and suppress 
the K dependence of the Jack polynomial by setting 

^ j'n\^,^)- (61) 

We now prove three statements 

Proposition 3 The asymptotic Bethe Ansatz equation JSutTWj 

N+l 

L I ' 2 



^ ( + ^ ^ sgn(fc, - %) I , heZ (62) 



for the quasimomenta ki are equivalent to the periodicity condition Prop. [I] of 
the creation (annihilation) operator a\ (au) ■ 

Proof In the region 

xi > X2 > ■ ■ ■ > xn+1 (63) 
the creation function a^(x, x') can be written as 

tr_M - /^nf<,Bin(.(x--x;)/L)\^ 



aT(x,x') = 



\Ylf<i^ sin(7r(a:i - Xj)/L)J 

(N N N N+l 

l[l[sin{n{x^^x'j)/L)l[ [] sin{n {x'^ - Xi}/L) 
j=ii=j+i 

(N+l N \ 

ikN+1 ^ Xj - ikN+i X! ■ (^'^) 



15 



After introducing complex variables 



f 2TTixi\ , f 2'Kix':\ , 

exp — ^ , z,' = exp — ^ (65) 



\ L 

Eq. (|27p becomes an integral representation for ^Ar+i(k, z) 
+ ^ ^ \/FTT \2m) Ji(N) ^ J 



(2z)-^^n^r'"'"^ n ^.^^^iv(k,z') , (66) 



with c?^[z'] = rii^i '^-^i- Using the form of Eq. (|55p for V'Af (k, z') as well as for 
ipN+i (k, z) on both sides of Eq. ((66)l 



. . _ g^W(2»)-^^ / i \^^^ £^^+(.+1)^ 



/^^^ d-[zV.(z, z') n (n, z') (67) 

•^^in i=l 

is obtained. We recall that /za(z,z') was defined in Eq. Q in the introductory 
section. Here n is the same partition of length N on both sides. Now periodicity 
°f 4iv+i(^'^') requires 

^N+i = TrkN+i + N{X + 1) , with un+i G ^ • (68) 

This equation can be iterated 

h = — (n,-(A + l)(i-l)) , (69) 

where we assume that all rii G N. We can restrict ourselves to positive in- 
tegers, since any negative integers can be absorbed by an appropriate boost. 
Subtracting the center of mass momentum K from every ki yields 

h ^ ^(^n, + ^^{N+l-2i)^ , l<i<N + l. (70) 

It is easy to verify that the ki in Eq. (jTO]) are solutions of the Bethe equation 
(l62|) . The integers li are related to n.i by li — ni + {N+l — 2i)/2. This completes 
the proof of Prop. [3] . 

Proposition 4 The action of the creation operator at is equivalent to the inte- 
gral representation for Jack polynomials, found by Olshanski and Okounkov in 
\0097l . 
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Proof In order to obtain the recursion formula for Jack polynomials as found 
in [0097] , we use Eq. (|68p and the well known property of Jack polynomials 



•^^("' ^) = ■^N^ ({"1 + 1: ^2 + 1, . . .}, z) = J^N\n + 1, : 



(71) 



Eq. ([57)1 is written as 



^^^^d^[z']MA(z,z')ji^\n-nAr+i,z') . (72) 

This is exactly the result by Okounkov and Olshanski f [0097| . Proposition 
6) for a Jack polynomial with partition n — rtAr+i, or equivalently for a Jack 
polynomial Jn.k-un+i (boosted by — tiat+i) with partition n, if we adjust the 
normalisation constants 

C^(k)(2z)-^^ — = [][B(n,,-nAr+i + (Ar+l-*)(A+l), A+l)]"! . (73) 

Here B{x,y) is Eulers beta-function. The quasimomenta ki on the 1. h. s. are 
related to the integers rii on the r. h. s. by Eq. ([M]) . This completes the proof 
of Prop, m 

Proposition 5 Let (n, z) be a Jack polynomial to partition n = {ni, . . . , n^} 
and let n be a positive integer. Define 

(74) 

T/ien i/ie following integral representation for Jack polynomials holds 

dz[ dz2... dz'^i^n,\{^, z) Jy {n,z') ^ BN6n,n,„JN^i{n~ n,z) . 

(75) 

For XeN, Um ^ ni + {X + 1){N -i), 1 <i < N. For X N, m ^ N. The 
contour integral is over the unit circle. The normalisation constant is 

N 

Bn = 2TTiY[[B{n,-n,n + {m~i){X + l),X+l)] . (76) 



i = l 



Proof The action of the annihilation operator Ok is defined by Eq. ([28|) . It is 
given in coordinate free notation by 



IV^(k)) = ^(-l)^+i4,fcjVjv-i(k^.)> , (77) 
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where in the state |'!/'7v_i(k^i)) a particle with quasimomentum ki has been 
deleted. After the variable transformation ((65|) using Eq. ([55]) yields 

/ d^[z>„,A(z,z')J^^^(n,z') , (78) 

where we have set 

n= ^k + {\ + l){N -1) . (79) 

ZTT 

The relation of the set of integers n to the set of quasimomenta k is given by 
Eq. ([69]) . Periodicity of afc(x, x') requires n G N. Therefore the Kronccker-delta 
for the quasimomenta transforms to 

SkM ~^ Sn,ni + {X+l){N-i) ■ (80) 

For i = N this is just (5„,„jv • For i ^ N it can only be non-zero for A G N. This is 
one assertion of Prop. [S] Due to the antisymmetry of the integrand in Eq. ((28l) . 
the lower bounds in I^^^ can be extended for all integration variables x[ to zero, 
respectively for all in Eq. (l78|) to one, without changing the integral. This 
yields Eq. (TTS]) and completes the proof of Prop. [H 



3.4 Rational Calogero—Moser— Sutherland system 

For the type (Il)-interaction potential the model is called rational CMS-model. 
Recursion formula (|27p was derived in a slightly different form in |GK02] . There, 
instead of ipN(k,x.), 

0Ar(k,x) = J] sgn(fc„ - fc^)sgn(a:„ - x,„) |^^(^^^^('|^)|A+i (^1) 

was considered. 0Ar(k, x) is completely symmetric in x and in k, as well as 
under interchange of the two sets x and k. For the special values A = —1/2, 0, 1 
it is the group integral ([6]) over U (TV), where U (N) is the unitary group over the 
real field (A = —1/2), the complex field (A = 0) or the quaternion field (A = 1). 
The parameter (3 of Ref. [GK02| is related to the coupling constant A by A = 
(3/2 — 1 (see footnote [Ij . The measure function 

/iA(x,x') - |A^.+l(x)r^-lat(x,x')|A^(x')|^+^ (82) 
has the geometrical interpretation as the invariant Haar measure over the coset 

OiN) ^ m 



N t 
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in a special parametrization, called Gelfand-Tzetlin coordinates [GT50[ IGK02] . 
Gi (g) G2 denotes the direct product group. Up to now there exist exphcit 
results for wave functions of the interacting rational CMS-model only for A = 1 
for small particle number up to = 4 and for three particles for arbitrary A 

pros] . 



4 Proofs 

In this section we prove Theorem [T] and Proposition [2] 

4.1 Proof of Theorem [T] 

We prove the five points of Theorem [1] 

1. In order to prove part [1] of Theorem [1] we show that the annihilation 
function afe(x, x') and the creation function a^(x, x') defined in (1^ are 
solutions of the differential equations ((2T|) . We first focus on the creation 
function a^(x, x'). We make for a[,(x, x') the Ansatz 



a]^(x, x') = cxp 



N+l 



N /N+l N \" 

n<m \n—l m—l / _ 



(84) 



where F{x) is an arbitrary even function. Acting with '^^^l -^rr and 
with J2n=i on a|.(x,x') yields 



^+1 Q2 



N 



N+l N 

)f{Xn ~ 2;;) - 2 E E ~ Xm)f(Xn - x'l) + 



1^711^71 

N N+l 

l^m n 
N+l N 



(85) 



4(x,x') 
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and, by the same token 



n=l ^-^ri 



N N iV+1 

7V+1 AT 

E E - ^!)/(^™ - -^J) - + (86) 

N+1 N 



4(x,x') 



where we defined f{x) = -^F{x) and /'(a;) — ^F{x). We now define 
the Hamiltonian in x representation, as 



N 



ij^(x) = -E 



92 



n— 1 n^m 

N 



(87) 



which comprises all terms in Eqs. (I85p and (j86|) . which depend only on 
one set of variables. Subtracting Eq. from yields 



i7Ar+i(x) -i7Ar(x')J al(x,x') = k^al{x,x') + 

N N+1 



JV+1 N 

E E [■^(^" ^ x[)f{Xm - x'l) + 2f(Xn - Xm)f{Xn ~ x'l)] 



4(x,x'). 



If f{x) fulfills the functional equation ([24|) 

f{x)f{y)+f{x)f{z) + f{y)f{z) = const. 

x + y + z = 0, (89) 

the term in the squared bracket on the left hand side becomes constani0. 



*With regard to Eq. I I88I I one might conclude that the functional equation 

/(xi - X2)f{xi -y) + f(x2 - xi)f{x2 -y) + f{xi - y)f{x2 ~ y) 
v{xi,X2) +u{y) 



(90) 
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We can rewrite Eq. ([55]) as 

HN+ii^) - i/Ar(x')] 4(x.x') = [e - 7V(iV + l)const.] 4(x,x') . 

(92) 

Observing that condition ([55)1 also yields 



^ 1 

J2 fixn-Xrn)f{xn-xi) = - -7V(7V - 1) (TV - 2)const. , (93) 



we find 



k^lix,^') = [HAr+i(x)-iJAr(x')]4(x,x') (94) 
^+1 Q2 



ri=l 

[/^(a;„ - Xm) - f'{xn - Xm) + const.] . (95) 

This is our assertion Eq. (|26p with dispersion relation 

e(fc) = fc2 . (96) 

The proof for the annihilation function afe(x, x') goes along the same lines. 
This completes the proof of part [T] and of part [2] of Theorem [T] 



2. In order to prove part [3] we have to show that the action of H^+i {Hn~i 
on the boundaries of the integration domain [lout ) does not give ris 
to boundary contributions. In other words, we have to show that 

i?7v+i(x 



iy^^^^d^[x>J(x,x')x^(x') = 

^^^^d^[x>S(x,x')i?^(x')xiv(x') 
HN-ii^) f d^[x>o(x,x')xiv(x') = 

d^[x']ao(x,x')i/jv(x')Xiv(x') , (97) 



with arbitrary functions v and u would be a sufficient, much weaker, condition on /. That this 
is not the case, is most easily seen by invoking translation invariance. Translation invariance 
of the left hand side of I I90I I requires translation invariance of the right hand side 

v{xi + a, X2 + a) + u(y + a) = v{xi,X2) + u{y) 

[— + — ]v{x^,X2) + —u{y) = 0. (91) 
V dx\ ax2 J ay 

Therefore we obtain v(x\,X2) = v{x\ — X2) and u(y) = const.. This means, the right hand 
side of Eq. I I90I I is independent of y. We now invoke permutation symmetry of the left hand 
side under interchanging x\ ^ y and X2 y- The same symmetry must hold on the left hand 
side, and therefore v{xi,X2) = const., too. Eq. I I90I I reduces to Eq. II89I I. 
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where we choose for convenience fcjv+i = 0. Here XAr(x') G 2? is a test 
function V cHn- The proof of Eq. ([W]) is given in App. S 

We now show that is an annihilation operator, i. e. 



afe|'(/'Ar(k)) = , iik^ki 



Cat 



(98) 



To this end (for the moment) we assume part[5]of Theorem[T]to be proven 
and consider the scalar product 



(7^^(k')|V^Ar(k)) 



rdet[27r(5(/c- - k.j)]i<i^j<N 



d[x]V'^(k',x)7/;Ar(k,x) 



d[x] 



d[x'] d[x"] [4, (x,x') 



V'^_i(k',x')4„(x,x")V'JV-i(k,x") 
c^[x']V'^~i(k', x') / d[x]afc/^ (x', x) 

^d[x"]4„(x,x")^iv-i(k,x") 

(V'jv-i(k')lafc^V^jv(k)) 
, if fc^ 7^ fci . . . /cat . 



(99) 



This proves Eq. 
I^iv-i(k')). 



since Eq. 



holds for an arbitrary wave function 



3. The pieces of part [4] of Theorem [T] which concern the action of the Hamil- 
tonian have been proven already before. To complete the proof we have 
to show in addition that a differential equation similar to Eq. (PTjl holds 
for the center of mass momentum operator PAr(x), defined in Eq. PT|) . 
namely 



[PAr+i(x)-P^(x')]4(x,x') = fc4(x,x') 

[^V(x) - Pjv+i(x')] afc(x,x') = -fcafc(x,x') . 
It is straightforward to see that this is true. 



(100) 



4. Finally we have to prove the orthogonality relation (i. e. part^of Theorem 
[T]). To this end we observe that the iV real numbers k = {fci,...,/cjv} are 
conserved quantities. This means that exactly TV mutually commuting 
selfadjoint operators In , 1 < n < N can be constructed with eigenvalues 
En — ^i- the potentials in Table [T] these operators can be 
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constructed with Dunkl operators Pm [Dun89[ IPol92j 



N 



m—1 



d 

-i^ 1- ^ f{Xm - Xn)P,ira , (101) 



where Pnm is the exchange operator defined in Eq. ([55]) . Therefore, for 
any two unequal sets k and k' at least one selfadjoint operator /„ can be 
found for which 

/„Vjv(k,x) ^ J„Viv(k',x) , ifk^k'. (102) 

Now we can invoke a fundamental theorem of functional analysis for self- 
adjoint operators: Two eigenfunctions of a selfadjoint operator to different 
eigenvalues are orthogonal. This completes the proof. 

4.2 Proof of Proposition [2] 

The five potentials of Table [1] have to be treated differently. 



Potentials (III)— (V) For the potentials (III)-(V) the normalisation constant 
CAr(k) is determined by the condition that in the asymptotic regime, where all 
particles are far apart, Xi — Xi+i oc M, M ^ oo, < N, the wave func- 
tion '(/)Ar+i(k, x) becomes the free Fermion wave function ■0^^^(k, x), defined 
in Eq. (|46)) . For the (5-interaction potential (V) the normalisation has already 
been calculated in Sec. 13.21 

In the asymptotic regime Af — > oo the cases of hyperbolic interaction poten- 
tials, (III) and (IV), can be reduced to the case of ^-interaction. To see this we 
introduce the rescaled particle coordinates Mxi and Mxi ~ x[ in the recursion 
formula Eq. ((27|) and obtain 

V'^-,i(k,Afx) = ^^^^"fl^ d^[x]4^^^(Mx,Mx)V'^(k, AfxI103) 

We now split the integration region into a region where all integration variables 
Xi are away from the boundaries Xj and the regions where one or more integra- 
tion variables are close such that (a;- — x^_^-^)/M ^ for M — ^ cxd. Since the 
integrand has by construction no singularities the second region has integration 
measure zero for Af — > oo and does not contribute in leading order in M . In 
the first region all integration variables Xi are away from the boundaries Xj and 
consequently apart from each other Xi — Xj = 0(M'^). Therefore we can substi- 
tute 'i/'Af (k, A/x) by the free particle wave function On the other hand the 
creation functions a[,^^^(Afx, Afx') for the three potentials (III)-(V) become 
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identical for M oo, since 

lim Aln| sinli(Maa;)| = lim A In | cosh(Af ax) | = Mc\x\ , c = aX. (104) 

A/— >oo M^oo 



« V .ys 

(III) (IV) 

In particular limM-.oo a\{M-x., Afx') = 1. We directly find, using Eqs. ([55)1 and 

^A,+i(k,Mx) = -V;^"^i(k,Afx) , (105) 

which is the result of Prop. [51 This completes the proof for potentials (III) to 
(V). 

Potential (I) and (II) For the trigonometric potential (I) the normalisation 
constant 



-N N 

\B{^{h^kN+l)A + l) 



C^(k)=fA) ^. (106) 



has been calculated already in Sec. 13.31 To obtain the normalisation constant for 
the rational potential (II) one can either follow the route of Sec. l4.2l and perform 
the integration in the asymptotic regime. But the evaluation of the resulting 
integral is by no means trivial. Therefore we resort to a different route. 

Eigenfunctions for the rational CMS Hamiltonian (II) emerge from the trigo- 
nometric case (I) in the limit L ^ oo keeping k and x finite. The same happens 
to their creation functions 

. ^ ^ ATA 

ai"'(x,x') = hm - a«(x,x'), (107) 

L^oo \ TT / 

where the upper index denotes the type of potential. This allows us to obtain 
the normalisation for the potential (II) case by taking the limit L ^ co oi 
Eq. p06p using the asymptotic expansion of the beta-function 

1 

lim — = (1 + lower order terms) . (108) 



X — ^oo 



The result is 

^^^^^ " r^(A + i) ni^'- W) ■ (109) 

It differs from the normalisation Eq. (5.7) of Ref. |GK02] by a /c independent 
factor. 
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5 Conclusions 



We constructed creation and annihilation operators for spinless interacting Fer- 
mions. Applying the creation operators successively onto the vacuum any N 
particle eigenstate can thereby be generated. The eigenstates are given as a 
{N — l)A^/2-fold integral. For the trigonometric CMS-Hamiltonian the equiva- 
lence of these eigenfunctions to other representations of the eigenstates has been 
demonstrated. 

The developed formalism paves a new way of searching and classifying inte- 
grable quantum systems, complementary or alternative to Bethe's Ansatz and to 
the Yang-Baxter equation. An interacting many-body Hamiltonian is exactly 
solvable if it can be transformed by a unitary transformation to a Hamiltonian 
containing one-body operators only. A method to use the present formalism 
for the calculation of thermodynamical quantities and of correlation functions 
is yet to be developed. 

It has to be stressed that the constructed operators have always Fermionic 
commutation relations. Therefore they differ fundamentally from the Bosonic 
operators which appear for instance in the Bosonisation approach. The latter 
describes the fundamental excitations of the N particle system by Bosonic op- 
erators which act on top of the filled Fermi sea. Thereby the Hamiltonian is 
effectively diagonalised. In contrast our diagonalisation is exact. 

The constructed operators are similar but not equal to the Fadeev-Zamolod- 
chikov operators [EK05J. The latter obey commutation relations which involve 
the two-body scattering matrix and need not necessarily have the commutation 
relations (|19p with the Hamiltonian. More on the connection of the creation 
(annihilation) operators, constructed here, to the Fadeev-Zamolodchikov alge- 
bra will be given in a subsequent publication [Koh| . where spin 1/2 Fermions 
are to be discussed. 

In the present work we focused on non-relativistic spinless Fermions. More- 
over translation invariance was assumed. It has to be stressed that the devel- 
oped formalism does not hinge on these assumptions. The extensions to spin 
1/2 Fermions and to Bosons will be given in separate publications. Applica- 
tion to lattice theories and 1-d relativistic field theories present an interesting 
challenge for the future. 

I acknowledge financial support by the German Research Council (DFG) 
with personal grant No. K03538/1-1 and within the SFB-TR/12. I thank 
F. Calogero, T. Guhr and A. Komnik for useful discussions. 

A Discussion of Boundary terms 

In this appendix some technical details concerning boundary contributions are 
reviewed. In particular the proof of Eq. (|97p is given in some more detail. 

Scattering Boundary Conditions (SBC) We look at the creation opera- 
tor, i. e. at the first line of Eq. ((97|) . The action of _ffAr4.i(x) onto the integral 



25 



yields after a somewhat cumbersome calculation 

J d^[x']iJAr+i(x)aJ(x,x')xjv(x') + boundary terms^^^ . (110) 

We considered for simplicity the operator aj. The results are easily extended 
to . The boundary terms are given by 



boundary terms 



,(1) 



5:/d^-Mx'^Jat(x',x)A^^(x' 
^ r f) 

5:yd^-Mx'^.]£-4(x',x)xA.(x' 

5:yd^-Mx'^J«S(x',x)£,x^(x' 

E yrf^-Mxy^«J(x',x)xAr(x' 



(111) 



We now use part 1 of Theorem [T] and substitute iJjv+i (x) by Hm (x') in 
Eq. piop . Integration by parts then yields 

j^^^^ d^[x']i/^(x')aS(x,x')XA^(x') - j^^^^ d^[x>t(x,x')ff^(x')xA'(x') 

+boundary terms*-^-* , (112) 

and altogether 

Hn+i{^) J^^^^ d^[x']aS(x,x')xiv(x') = 

y^^^^^d^[x>t(x,x')i/^(x')xiv(x') 

+boundary terms^^-* + boundary terms'-^-' . (113) 



This is the assertion of Eq. (|28p if the two boundary terms cancel each other. 
A direct calculation, using the identity 



d 



= -^4(x',x) 



(114) 



shows that this is so. 
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We now look at the action of the annihilation operator, i. e. at the second 
line in Eq. (jOT]) . Acting with i/w-i onto the integral yields an expression similar 
to Eq. (flT^ 

i?jv-i(x) / d^[x']ao(x,x')xiv(x') = 

d^[x>o(x,x')ffjv(x')Xiv(x') (115) 

(N) 
out 

+boundary terms'^'' + boundary terms'^^^ . 

However, now the two boundary contributions do not cancel each other iden- 
tically. The boundary contributions at ±00 vanish if the annihilation function 
ao(x, x') is bound from above such that for xn{^') G D the product ao(x, x') 
XAr(x') is integrable: 

3C e R : ao(x, x') < C , for x'l — > 00 and x'j^ — > — cx) . (116) 

For all cases, considered in Table [TJ this is indeed the cas^. 



Periodic Boundary Conditions (PBC) For the action of the creation op- 
erator the treatment of the boundary contributions is completely analogous to 
the treatment for SBC. Here, in addition, we observe that periodicity of the 
integral 



/Ar+i(fc,x) = 




(117) 



in all unprimed variables Xi, 1 < i < N + 1 with period L requires on the one 
hand a^(x, x') to be periodic in x.^, 1 < i < + 1 and on the other hand 

L 

dx^4(x,x')xjv(x') = , yi<i<N. (118) 

In order to see that for the action of the annihilation operator no boundary terms 
occur, we rewrite the integration domain I^^^ in the second line of Eq. (j97p as 



/Af_l(fc,x) = 




= / dx[ I dx2... da;^afe(x,x')xAr(x') . (119) 

Jo Jo Jo 

The second equation was obtained due to the total antisymmetry of the inte- 
grand. Now it is seen that periodicity of /Ar_i(x) in all variables Xi, 1 < « < 

^For particles with (5-intcraction this is only so, because we are considering spinless fermions 
and ao{x, x') = 1, see Sec. 13.21 
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N —1, requires the integral over one period to vanish for all integrations but for 
the integration over x'l 

L 

dx^afc(x,x')xiv(x') = , Vz^l. (120) 
This allows us to write I]y^i{k,x.) as 

/Ar_i(fc,x) = / dx'i I dx'^... / dx'j^ TT 

isgni(a;, -x^+i)afc(x,x')xA'(x') , (121) 

where sgn^(a;) is the periodic sign-function sgn^(x) = sgn^(a; + L). In this form 
it is clear that 

-ffAr-i(x)/Ar_i(fc,x) = / dx[ dx'^... dx'pf Yi o^S^Li^i - ^t+l) 

Jo Jo Jo ^ 

afe(x,x')ffjv(x')xiv(x') (122) 

if the integrand is a periodic function and in particular if afc(x, x') is periodic 
with period i in all x^, 1 < i < A^. This is the assertion of Prop. [TJ Indeed, 
for PBC the form (I12ip serves as a definition of the integration domain of the 
annihilation operator alternative to ij^^^ defined in Eq. 
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